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Abstract 

Twist grain boundaries are widely observed in lamellar phases of block copolymers. A mesoscopic 
model of the copolymer is used to obtain stationary configurations that include a twist grain 
boundary, and to analyze their stability against long wavelength perturbations. The analysis 
presented is valid in the weak segregation regime, and includes direct numerical solution of the 
governing equations as well as a multiple scale analysis. We find that a twist boundary configuration 
with arbitrary misorientation angle can be well described by two modes, and obtain the equations 
for their slowly varying amplitudes. The width of the boundary region is seen to scale as e -1 / 4 , 
with e being the dimensionless distance to the order-disorder transition. We finally present the 
results of the linear stability analysis of the planar boundary. 

PACS numbers: 47.54.-r.61.25.H-. 83. 50.-v.05.45.-a 
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I. INTRODUCTION 



Block copolymers are being explored for either direct use or as processing templates in 
nanolithography, photonic devices, high density storage systems, or drug delivery, just to 
name a few examples [l, 2, 3, 4]. Self assembly into mesophases of different symmetries and 
of controllable periodicity (typically at the nanoscale) makes these materials a very versatile 
tool, and hence the interest in studying their architectures and self assembly mechanisms 
[l|, Q]. A common practical limitation to widespread use, however, is the considerable dif- 
ficulty encountered in producing well ordered microstructures jsl, [f], ?]]. Given that the 
longest relaxation times of partially ordered microstructures are often controlled by existing 



topological defects, much attention has been paid to the motion of disclinations 
grain boundaries in lamellar jl], n]] and cylindrical phases 11]. 



and 



Theoretical analyses of defect motion have been based on a mesoscopic description of a 
copolymer melt which is valid for characteristic time scales much longer than the slowest 
relaxation time of the polymer chain jl2,[l^, Q]- Asymptotic methods commonly employed 



in studies of defect dynamics in systems outside of equilibrium [15( have been applied, to tilt 
grain boundaries in lamellar phases Q, Q • This type of boundary separates two domains 
of differently oriented lamellae such that the plane formed by lamellar normals of the two 
domains is perpendicular to the boundary plane. Examples include boundary mig ration 



induced by lamellar curvature {9], Q, [iSl] and the effect of an imposed shear flow |20|, [21]. 

In three-dimensional samples, 90° tilt boundary configurations (the so called T-junctions) 
are rarely observed in experiments 22|, possibly because they are generically unstable 23]. 
On the other hand, twist boundaries (such that the wavevectors of both adjacent lamellar 
domains lie on the boundary plane (see Fig. [T|)) of various misorientation angles are com- 



monly observed 



very limited 



24 
28 



25 . [26|. Nevertheless, analyses of their structure and stability are still 
29]. 



We focus in this paper on a coarse grained model of a twist grain boundary, leading 
to the associated amplitude (or envelope) equation description. We obtain a stationary 
profile comprising a twist grain boundary, and numerically compute its linear stability. Our 
results are based on the Leibler or Swift-Hohenberg model [l^, S]> valid in the limit of weak 
segregation. The analysis is conducted for a boundary of arbitrary misorientation angle a. 
In contrast with the results obtained for the case of tilt grain boundaries, we find that the 
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FIG. 1: Schematic of a twist grain boundary separating two lamellar domains with a misorientation 
angle a. Also shown (bottom) is the morphology at the boundary as given by the order parameter 
of the model. 

twist boundary width scales as e~ 1//4 , with e the dimensionless distance to the order- disorder 
threshold, and that the twist boundary is linearly stable to long wavelength modulations for 
any angle a, consistent with experimental findings in copolymer melts. 



II. MODEL 



A. Coarse-grained model equation 



At a mesoscopic level, a weakly segregated diblock copolymer melt close to the order- 
disorder transition temperature Tqdt is described by a free energy, function of monomer 



composition, given by Leibler 
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31 



the Swift-Hohenberg model equation 



14 



he corresponding relaxational dynamics leads to 



15 



30] . For a symmetric diblock melt (with equal 



volume fraction of the two constituent monomers), this model equation is (in dimensionless 
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units) 

^ = 6^ - (V 2 + ql) V - V 3 , (1) 

where the order parameter field ip represents the local density difference between the two 
monomers of the diblock, and go = 1 after rescaling, although we retain the symbol go i n 
what follows for clarity of presentation. As already stated above, the control parameter e 
measures the distance from the order-disorder transition or bifurcation point at which e = 0. 
For e > (temperature below Tqdt), a pattern with lamellar symmetry emerges, although 
usually accompanied with large amount of defects. Therefore typical configurations display 
a multidomain microstructure. 



B. Amplitude equations of twist grain boundaries 



Following the standard multiple scale approach in the weak segregation limit [15|, Il7| . 
we can separate the fast spatial/temporal scales of a base lamellar pattern from its slowly 
varying amplitude, and derive the associated amplitude equations for a twist grain boundary. 



The derivation is based on the method given in Ref. 
expanded as the superposition of two base modes 

1 



2l|. The order parameter field ip is 



1> 



[A exp (iqi ■ r) + B exp (iq 2 ■ r) + c.c. 



(2) 



where q\ = q x and q 2 = qo (cos ax + sin a y) (with a the twist angle) are the orientations 
of two domains adjacent the twist boundary (see Fig. [1]). The evolution of the complex 
amplitudes A and B is governed by (to leading order in 0(e 3 ^ 2 )) 



d t A 
d t B 



e-(Vl + 2iq d Hi y 
e- (Vl +2iq d ri2 y 
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-2\A\ 2 B 



(3) 
(4) 



where n^n^ are the normals to the lamellar planes in domains A and B respectively, 

is the the Laplacian operator on the lamellar plane of domain A, and V| 2 represents the 

Laplacian operator on the lamellar plane of domain B. For instance, if n\ — x (i.e. q\ = qox), 



d 2 + d 2 and d Hl = d x . 
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III. STATIONARY TWIST GRAIN BOUNDARY CONFIGURATION 



The steady state configuration of twist grain boundaries has been first examined through 
direct numerical solution of the model equation ([T]). A pseudospectral method in Fourier 
space is adopted, with periodic boundary conditions in all three directions. A Crank- 
Nicholson time stepping scheme is applied to the linear terms, with a second order Adams- 
Bashford algorithm used for the nonlinear term. Periodic boundary conditions are satisfied 
through the consideration of an initial configuration comprising a symmetric pair of twist 
boundaries that are sufficiently far apart so that their motion is approximately independent. 
An additional restriction needs to be placed on the dimension of the computational cell along 
the x and y directions on the plane of the grain boundary (as shown schematically in Fig. 
[2]). Due to the requirement that an integer multiple of lamellar periods Ao (= 2ft /qo) must 
equal the length of the computational cell in the direction parallel to the lamellar normal, 
the unit cell lengths are l x = A / sin(a/2) and l y = A / cos(a/2). We consider a uniform 
spatial discretization L z grid with spacings Ax = Ay = l y /16, and 

Az = Ao/16, corresponding to 16 grid points per unit cell length. Most calculations shown 
below correspond to a system size of 256 3 , with a dimensionless time step At = 0.2 used in 
the numerical integration. 

A typical stationary configuration is shown in Fig. [31 corresponding to a = 90° with 
e = 0.04 and at a time t = 10 4 . Also shown in the figure (in grey scale) is the two 
dimensional order parameter at the boundary interface. It is doubly periodic along the two 
directions defined by the bulk lamellar domains adjacent to the grain boundary. 

Before we carry out a multiple scale analysis, we have checked the underlying assumption 
that the stationary order parameter field ip can be decomposed into two Fourier modes (see 
Eq. ([2])). We calculate the Fourier spectrum of the order parameter field both at the grain 
boundary and in the bulk. We illustrate our findings with the case of a a = 75° twist 
boundary with e = 0.02. The two dimensional Fourier spectrum of ip on the boundary plane 



(z = (59/256)L 2 ) shows four maxima at wavevectors (±q X Q,±q y o), with y ql + q y0 = qo 
(q is the wavenumber in the bulk). We also observe 2 arctan(g a , / ( ?yo) — 75° (exactly the 
misorientation angle a). This reflects the fact that the order parameter in the grain boundary 
region is a combination of the two bulk modes. The same conclusion is supported by an 
analysis of higher harmonics in the spectrum. Figures and Hb show the intensity of the 
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FIG. 2: Schematic of the constraint imposed by periodic boundary conditions. Solid and dashed 
lines represent two misoriented lamellar layers at the grain boundary. In order to accommodate 
two domains with a relative misorientation a and use periodic boundary conditions, one needs to 
choose as unit cells dimensions l x = Ao/ sin(a/2), l y = Ao/ cos(a/2) along the x and y directions 
on the boundary plane. 

spectrum along q x at q y = q y o and at two different values of z: one within the grain boundary 
(Fig. Hk), the other in the bulk A phase (Fig. Hb) (identical conclusions can be drawn from 
the analysis in phase B). Figures Hb andHJi show the same quantity but as a function of q y 
at q x = q X Q, and for the same two values of z. The fact that all the visible harmonics within 
the grain boundary region are almost the same as in the bulk suggests that in the weak 
segregation limit considered here, the superposition of the two bulk modes in Eq. ([2]) used 
for the multiple scale analysis that will follow appears to be sufficient for the description of 
the order parameter profile around the grain boundary. 

IV. ASYMPTOTIC BOUNDARY WIDTH 

The boundary width 5 of a twist grain boundary as a function of e in the limit e — > 
can be determined either numerically from the stationary configuration given above, or via 
a multiple scale analysis ([3]) and (j4j). In the latter case, simple dimensional analysis of Eqs. 
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FIG. 3: A stationary a = 90° twist grain boundary configuration in a lamellar phase, as given by 
numerical solution of the model equation ([I]). The grid size is 256 3 , and e = 0.04. Right panel: 
order parameter (gray scale) at the boundary interface. 



(a) z = 59/256L 



(b) z = 5/256L z 




intensity 
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FIG. 4: Power spectrum of ifi along specific directions in q = (q x , q y ) space, and at constant location 
z for a a = 75° twist boundary. We show the spectrum for e = 0.02 and t = 2000. Panel (a) shows 
the power spectrum as a function of q x for q y = q y o and z = (59/256)L 2 , whereas (b) shows the 
spectrum far into the bulk at z = (5/256)L z . Panels (c) and (d) show similar spectra as a function 
of q y for q x = q X Q and at z = (59/256)L z (c) or (5/256)L 2 (d). Here (q x o,Qyo) is the location of the 
peak of the two dimensional power spectrum for the order parameter tp. 
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and (jlj) along the grain boundary normal (the z direction) leads to the following result 

5 ~ e- 1/4 . (5) 



This is in contrast with the known behavior for a tilt gain boundary in which 5 ~ e -1 / 2 [16]. 
The latter scaling behavior follows from the fact that there are two distinct characteristic 
length scales for lamellar relaxation: One along the direction parallel to the lamellar normal 
with scale l± oc e -1 / 2 , the other along the plane of the lamella with scale l\\ oc e -1 / 4 . For a 
twist grain boundary, on the other hand, the direction (z) normal to the boundary is parallel 
to the lamellar planes of both phases, and hence it is reasonable to expect that the boundary 
width along z scales as e" 1 / 4 . In summary, a twist boundary is much narrower than a tilt 
boundary in the limit e — > 0. 

The result above holds for any misorientation angle as we have verified by numerical 
solution of the model equation ([1]). We first determine the location of the boundary by 
estimating the amplitude B(z) |2l| 

B{z) ~ — ^2 ftH^" ™b = m\ ; z) - ip(r- h B = (m - 1/2)A ; z)\ , (6) 

m=l 

with fiB the unit normal to lamellae B and N the number of pairs of lamellae. The boundary 
region is chosen such that the value of B(z) lies within 10% - 90% of its maximum. Since 
the width of the boundary is only several times the lamellar width, a linear interpolation 
algorithm is used to increase the accuracy of the boundary location. The relation obtained 
between boundary width 5 (in dimensionless units) and misorientation angle a is plotted in 
Fig. [5] for e = 0.02. For a > 20°, the boundary width becomes approximately independent 
of a. Otherwise, 5 increases rapidly with decreasing angle. Although the accuracy of our 
numerical solution degrades when a is small, the trend obtained points to a divergence 
of the boundary width as a — > 0. We find similar results when numerically solving the 
corresponding amplitude equations and (j4j) for small twist angles. Figure [5] also shows 
our results for 5 as a function of e for a = 90°. Given the spatial discretization used in our 
integration, the boundary widths that we have been able to investigate range from 3Ao at 
e = 0.001 to 3Ao/4 at e = 0.4 (with Ao = 2n/q ). Within this limited range, a power law 
dependence between 5 and e is found, with an exponent —0.244 ± 0.002, in agreement with 
our expectation from dimensional analysis. Analogous results have been obtained for other 
values of a. 
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FIG. 5: Boundary width 5 (in dimensionless units) as a function of twist angle a for e = 0.02. 
Inset: Width 5 vs. e as obtained from the stationary solutions of Eq. (JTJ) for a = 90°. The slope 
of the log-log plot is -0.244 ± 0.002. 



V. STABILITY ANALYSIS 



As noted above, twist grain boundaries are observed in great abundance in 



experiments 



that address the microstructure of large samples in the lamellar phase [24j, [25|, [26|]. We 
conduct here a linear stability analysis of a planar boundary from the amplitude equations 
([3]) and (jl]) that are derived from our model equation. We start from a base state involving 
a stationary and planar twist boundary of arbitrary misorientation angle a and wavenumber 
Qo |32j- The corresponding amplitudes and are assumed to be only a function of 
z, the direction normal to the boundary, and are given by 

eA® - %A<® - \A®\*A® - 2\B^\ 2 A^ = 0, (7) 
eB (o) _ a* B (o) _ |£(o)| 2jB (o) _ 2 \A^\ 2 B^ = 0. (8) 
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We next expand the complex amplitudes around the stationary solutions 

A(x,y,z,t) = A®(z) + J2Mq*,Qy,z,t)e i ^ + ™\ (9) 

B{x,y,z,t) = B^(z) + J2B(q x ,Qy,^t)e^ + ^\ (10) 

qx,q y 

substitute these expansions into Eqs. ([3]) and (jlj), and linearize the resulting equations with 
respect to the perturbations A and B. We find 

d t A(q x , q y , z, t) = [e - (d 2 z - q 2 y - 2q q x ) 2 - 2\A^\ 2 - 2\B^\ 2 } A(q x , q y , z, t) 

- (A^) 2 A*(-q x , -q y , z, t) - 2A^B^*B(q x , q y , z, t) 
-2A^B^B*{-q x ,-q y ,z,t), (11) 

d t B(q x ,q y ,z,t) = [e - (d 2 - q 2 y2 - 2q q X2 ) 2 - 2|^°)| 2 - 2\B™\ 2 } B(q x ,q y ,z,t) 

- (£?W) 2 B*(-q x , -q y , z, t) - 2B^A^*A(q x , q y , z, t) 
-2A^B^A*(-q x ,-q y ,z,t), (12) 

where q X2 = cos a q x + sin a q y and q y2 = — sin a q x + cos a q y . 

Since we do not have an analytic expression for the amplitudes of the base state, we study 
its stability by examining the temporal evolution of small random perturbations to both real 
and imaginary parts of A and B for a range of values of q x and q y , and integrating the system 
of Eqs. (!71)- (fT2l) numerically. The details of the numerical algorithm and procedure are given 
in Ref. [23|. The parameters chosen here are Az = Ao/8 for the discretization along the 
z direction, with L z = 1024 grid nodes (or equivalently a length of the computational 
domain of 128Ao). The time step chosen is At = 0.2. If the planar grain boundary is 
stable, perturbations in A and B will decay in time for all wavevectors (q x , q y ); otherwise an 
instability would manifest itself by an increase of these perturbations within a certain range 
of wavevectors. 

From the relaxation of the perturbations, we estimate the perturbation growth rate 
&(q x ,q y ) from |-B(t)| oc e a ^ 9x,9y ^ t and the numerical solutions for A and B for a given 

set of {q x ,q y ). A typical result is shown in Fig. E] for a = 90° and e = 0.04. We always 
observe that a < 0, for all the wavevectors of the perturbation explored. This is also the 
case for different values of e and angle a, as shown in Fig. [71 The maxima of a for a ranging 
from 30° to 90°, and e from 0.005 to 0.08 have been calculated, all yielding a stable planar 
boundary. 
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FIG. 6: Perturbation growth rate a as a function of wavevector (q x , q y ) for e = 0.04 and a = 90°. 
We find that a < over the whole range of wavevectors investigated. 




FIG. 7: Maximum perturbation growth rate <r max as a function of e, for different twist angles 
a = 30°, 45°, 60°, and 90°. 

VI. CONCLUSIONS 

We have used the Swift Hohenberg model as an approximate mesoscale description for 
the evolution of a twist grain boundary in a lamellar phase of a diblock copolymer. We 
have shown that the order parameter field can be well approximated in the weak segregation 
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regime by a combination of the two modes of the ordered lamellar phases on either side of the 
grain boundary. The equations governing the slow evolution of the amplitudes or envelopes of 
these modes have been derived for arbitrary misorientation angle. The stationary solution is 
only a function of the coordinate normal to the grain boundary plane, and is characterized 
by a width 5 ~ e~ 1//4 , with e the distance from the order- disorder point. We have then 
conducted a linear stability analysis by direct numerical solution of the governing equations, 
and have found that the twist boundary is linearly stable within a wide range of parameters 
investigated. 
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